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Abstract 

We explain a profound complexity of magnetic interactions of some tecimologically relevant 
gadolinium intermetallics using an ab-initio electronic structure theory which includes disordered 
local moments and strong /-electron correlations. The theory correctly finds GdZn and GdCd 
to be simple ferromagnets and predicts a remarkably large increase of Gurie temperature with 
pressure of -1-1.5 K kbar“^ for GdGd confirmed by our experimental measurements of -I-1.6 K 
kbar“^. Moreover we find the origin of a ferromagnetic-antiferromagnetic competition in GdMg 
manifested by non-collinear, canted magnetic order at low temperatures. Replacing 35% of the 
Mg atoms with Zn removes this transition in excellent agreement with longstanding experimental 
data. 
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Lanthanide compounds play an increasingly important role in the development of novel 
materials for high-tech applications which range from mobile phones and radiation detec¬ 
tors to air conditioning and renewable energies. Much of this stems from their magnetic 
properties, so that they are indispensable components in permanent magnets, [1] magneto- 
responsive devices for solid state cooling [2] and other applications. Common to all the 
lanthanide elements is their valence electronic structure which makes them chemically sim¬ 
ilar and also causes magnetic order. Lanthanide atoms are predominantly divalent (5d°6s^ 
valence electron conhguration), becoming mostly trivalent in a solid, donating three va¬ 
lence electrons to the electron glue in which the atomically-localised /-electron magnetic 
moments sit. The interaction between these moments derives from how the electron glue is 
spin-polarised. The longstanding RKKY [3] (Ruderman-Kittel-Kasuya-Yosida) free electron 
model of this electronic structure is typically used to try to explain the many features of 
the indirect coupling of the 4/-electron moments despite its rather poor representation of 
the narrow band Sd-states. The possible importance of the latter has already been inferred 
from some earlier electronic structure studies gHT]. 

Whilst theoretical aspects of lanthanide magnetism are well understood at the phe¬ 
nomenological level, predictive hrst principles calculations are challenging owing to the com¬ 
plexities of the strongly correlated /-electrons and itinerant valence electrons along with the 
magnetic fluctuations generated at hnite temperatures. In this letter we explore lanthanide 
compounds with an ab-initio theory based on Spin Density Functional Theory (SDFT) in 
which the self-interaction corrected (SIC) local spin density (LSD) method [SI H] provides 
an adequate description of f-electron correlations unHn] and the disordered local moment 
(DLM) theory [13] handles the magnetic fluctuations. We are able to give a quantitatively 
accurate description of the diverse magnetism of Cs-Cl (B2) ordered phases of Gd with Zn, 
Cd and Mg which we test against experimental data and show the complex role played by 
the spin-polarised valence electrons. 

Local moments of hxed magnitudes are assumed to persist to high temperatures and 
in lanthanide compounds are formed naturally from partially occupied localised 4f-electron 
states. The orientations of these moments fluctuate slowly compared to the dynamics of 
the valence electrons glue surrounding them. By labelling these transverse modes by local 
spin polarisation axes hxed to each lanthanide atom i, bj, and using a generalisation of 
SDFT [i3](+sic m ng) for prescribed orientational arrangements, {bj}, we can deter- 


2 


mine the ab-initio energy for each conhguration, [T5HT9] so that the conhguration’s 

probability at a temperature T can be found. The magnetic state of the system is set 
by an average over all such conhgurations, appropriately weighted, and specihes the mag¬ 
netic order parameters, {nij = (bj)}, where the magnitudes rrii = |mj| range from 0 for 
the high temperature paramagnetic (PM) (fully disordered) state to 1 when the magnetic 
order is complete at T = OK. A distribution where the order parameters are the same 
on every site, {nij = m^z} say, describes a ferromagnetically ordered (FM) state whereas 
one where the nij alternate layer by layer between iriaX and —niaX characterises an anti¬ 
ferromagnetic (AFl) order. The free energy function J^({mj}), written in terms of these 
magnetic order parameters, nij, monitors magnetic phase transitions. It contains the effects 
of the spin-polarised valence electronic structure which adapts to the type and extent of 
magnetic order [T8l 12011^ . For lanthanide materials DLM theory describes how valence 
electrons mediate the interactions between the f-electron moments. These can turn out to 
be RKKY-like, but can also show strong deviations from this picture as we hnd here for 
simple Gd-containing prototypes. 

We start with GdZn, of particular interest in solid state cooling, ra but also because we 
expect its electronic structure to be straightforward PEa]. The Gd atoms occupy a simple 
cubic lattice of the GsGl(B2) ordered phase. Our hrst-principles SIG-LSD calculations hnd 
the ground state Gd-ion conhguration to be trivalent (Gd^^), with seven localised f-states 
constituting a stable half-hlled shell, in line with Hund’s Rules [TUI im [TU] . So Gd of all the 
heavy lanthanides has the relative simplicity of an S-state for its f-electrons, largely uncom¬ 
plicated by crystal held ehects and spin-orbit coupling. This permits a clinical look at how 
the interactions between large 4f-magnetic moments are mediated by the valence electrons. 
These come from both the lanthanide (5d^6s^) and the post-transition metal Zn which has 
low-lying, nominally hlled, d-shells (3d^°) added to its two s-electrons. Our ab-initio DLM 
theory can thus investigate the ehect of the lanthanide 5d electrons hybridising weakly with 
3d states. This touches on a very important aspect of many magnetic materials contain¬ 
ing both rare earth and transition metal elements PI where understanding the interplay 
between the localised lanthanide magnetic moments and the more itinerant magnetism orig¬ 
inating from the transition metal d-electrons is paramount for the design of more efficient 
materials. 

Our DLM theory calculations for the paramagnetic state of GdZn produce local moments 
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of magnitude /i ~ 7.3/iB on the Gd sites pointing in random directions so that there is no 
long range magnetic order, {rrij = 0}. The calculated paramagnetic susceptibility [151 dSl 
[TS] , x(q)) with a maximum at wave-vector q^ax = (0,0,0), shows that, in accord with 
experiment, GdZn develops ferromagnetic (FM) order below a Gurie temperature, Tc = 184 
K (at theoretically determined lattice constant, a = ath = 6.62 a.u.), somewhat lower than 
the experimental value of Tc = 270 K [25] (at a = Qexp = 6.81 a.u.). We hnd that GdZn’s 
Tc gradually decreases under pressure, P, with calculated = -0.45 K kbar“^ which 
agrees reasonably well with experimental value of -0.13 K kbar“^ from the literature [ 25 ] 
(Fig. 1(a)). The negative is typical of many metallic magnets owing to pressure-induced 
band broadening and diminished energy beneht from spin polarising the valence electrons 
around the Fermi energy. 

Naively one might expect similar effects if Zn is replaced with isoelectronic Gd whose hlled 
4(i—band states are simply more extended than the 3(i’s of Zn. Our calculations, however, 
show something rather different. Whilst both theory and experiment hnd GdGd to be a 
simple ferromagnet like GdZn, with Tc = 234 K {ath = 6.98 a.u.) and 265 K[2Z](aea;p = 
7.09 a.u. [IE]), in sharp contrast to its results for GdZn, theory predicts its Tc to increase 
quite dramatically with pressure (Fig.l), i.e. a positive and rather large Owing to the 
paucity of reliable published experimental pressure data for GdGd, [28] we have carried out 
measurements [18] to test this specihc prediction and a comparison between the calculated 
and experimentally observed Tc’s for GdGd as a function of pressure is shown in Fig.1(a). 

The theory-experiment agreement is excellent: ^ from theory is -1-1.5 and from exper¬ 
iment is -1-1.6 K kbar“^. Whilst not unusual for hrst-order magnetostructural transitions 
(e.g. ~ 1-3 K kbar“^ is observed in Gd 5 Sia;Ge 4 _a; alloys IZSl), this is a rather high rate 
for a second-order transition as occurs in GdGd. Reasons for this stark difference between 
GdZn and GdGd are found from our Tc calculations as a function of lattice parameter, a, 
(Fig. 1(b)). Starting from large values, Tc initially increases with decreasing Gd-Gd distance 
for both GdZn and GdGd, reaching a maximum whence it starts decreasing with further 
reduction of the Gd-Gd distance. The shown in Fig. 1(a) originate from where the two 
compounds have their equilibrium lattice spacings, ath, marked by blue (GdZn) and red 
(GdGd) arrows in Fig. 1(b). 

For materials with the same number of valence electrons per atom, the RKKY account of 
magnetic interactions would be the same. GdMg is isoelectronic with both GdZn and GdGd 
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FIG. 1: (a) Comparison between theory (open symbols) and experimental [18] (filled symbols) Tc 
differences, (T’c(P) — Tc(0)), as a function of pressure, P, for GdZn (blue circles) and GdCd (red 
squares). The experimental data for GdZn are from Ref. (|26|) . (b) Tc of GdZn (blue circles) and 
GdCd (red squares) as a function of lattice parameter a (atomic units) calculated from the theory. 
The vertical arrows indicate ath, red for GdCd and blue for GdZn. 


but with no filled 3d or 4d band of states. This difference leads to our DLM theory hnding 
GdMg’s PM state to be qualitatively different than GdZn’s and GdGd’s. We hnd a discor¬ 
dant blend of FM and AFl dominant magnetic correlations in the PM state - the calculated 
paramagnetic x(q) has two comparable peaks at wavevectors (0,0,0) and (0,0, |) [TH] (in 
units of ^). Which one is stronger depends on the a values used. At the theory volume 
{ath =7.00 a.u. [c.f. Qexp = 7.20 a.u. [30|]), our calculations predict a FM state below 
Tc =128K. Reducing the Gd-Gd separation weakens the FM aspects and, for example, a 4% 
decrease leads to an AFl state instead, below the Neel temperature T/v =87K. 

We determine the magnetic order that evolves as T is lowered through the transition 
temperatnre to OK as a conseqnence of these competing FM and AFl effects by using onr 
DLM theory [18] for the hrst time to describe a magnetically ordered state with a canted 
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structure and repeating the analysis for a number of a values. We set the order parameters, 
riij’s for the system at various stages of partial onto complete magnetic order, to alternate 
between mfZ + niaX and nifZ — niaX on consecutive Gd layers along the (1,0,0) direction 
giving a canting angle between layers of ©c = 2 arctan(^) so that the overall magnetization 
of a system is local moment size /i (7.3 /xs) times rrif. rrif ^ 0, rua = 0, ©c = 0 signihes a 
FM state and rUa 7 ^ 0, m/ = 0, ©c =180° a AFl state. 



FIG. 2; (a) The magnetic phase diagram of GdMg, represented by the canting angle Qc{T) and 
its dependence on lattice spacing alongside a schematic picture of the CM state, (b) The magnetic 
phase diagram of GdMg(i_ 3 ,)Zn 3 ;, 0c(T), and its dependence on x for a fixed lattice spacing equal 
to the 2% reduction value in (a), equal to ath of GdMgo, 6 Zno, 4 . 


Fig.2(a) summarises our results. These are the hrst ab-initio calculations to show canted 
magnetism (CM) in GdMg. The hgure shows the emergence of a CM from either a FM 
(©c = 0) or AFl state (©c =180°). For low T, ©c ranges from 70° at the theoretical 
equilibrium lattice constant (0% reduction) through to 120° (4% reduction) before eventually 
forming an AFl magnetic structure (angle 180°) with further reduction. This agrees very 
well with experiment [3T] which hnds that, upon lowering the temperature, GdMg orders 
into a FM state at ~ IlOK and then undergoes a further second order transition into a 
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canted magnetic ordered state at Tp 85K [3T1[32]. At low T the magnetisation 5 /ip, a 
valne we have also conhrmed with onr own experimental measnrements. This is indicative 
of the FM and AF components, ruf and nia, being roughly the same size giving a canting 
angle between 7 /ip-sized Gd moments of roughly 90°. This state is robust against applied 
magnetic helds [32| of up to 150 kOe. The experimental results are matched almost exactly 
by our calculations shown in Fig.2(a) for a 2 % lattice spacing reduction from ath- Liu et 
al. [31] also found that under pressure GdMg orders into a AFl from a PM phase at ~ lOOK 
and at a lower temperature undergoes a further hrst order metamagnetic transition into a 
canted FM phase. The authors estimated the pressure derivative of the magnetisation to be 
-0.04 fiB kbar“^ at 4.2K, which we have also conhrmed experimentally and in fair agreement 
with our calculated low T value of -0.02 /xp kbar“^. 

Experimentally it is known that when Gd is replaced by Tb in GdMg, there is a 1% 
lattice contraction [39| and a FM state undergoes a transition into a canted magnetic struc¬ 
ture at low T with 0^ of at least 90°. Replacing Gd with Dy leads to a larger lanthanide 
contraction and measurements jlOl show that DyMg orders into an AFl state, developing 
non-collinear structure with a FM component at low T and ©c of about 110°. This cor¬ 
relates with Fig.2(a) [H] for the smaller lattice spacing regime. The little available data 
for Ho-Mg [39] also indicates canted AF magnetic structure at low T. So we infer that 
the lanthanide contraction [15] in part causes the transition from FM-canted to AF-canted 
magnetic structures as the heavy lanthanide series is traversed. Our hgure 2(a) also implies 
a tricritical point (PM-AF-FM) at some concentration, y, in the (Tbi_yDyj^)-Mg alloy sys¬ 
tem with a transition to a canted structure at a marginally lower temperature or possibly a 
transition into a canted structure directly. 

This unusual canted magnetism of GdMg is evidently destroyed by nominally hlled, low- 
lying 3d or 4d bands from the non-lanthanide constituent. Our calculations, Fig.2(b), show 
what happens when a fraction x of the Mg sites in GdMg is replaced by Zn. T^. increases 
with X, and the low temperature canted structure vanishes altogether for x > 0.35. This ob¬ 
servation is in excellent agreement with the experimental data for GdMg(i_a,)Zna; of Buschow 
et al. [33] who gave an early report of a serious shortcoming of the RKKY picture. 

The successful capturing of these unusual temperature and pressure trends of the Gd 
intermetallics’ magnetism is a consequence of the theory’s detailed description of the valence 
electrons. The theory includes both the response of these electrons to the magnetic ordering 
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FIG. 3: The local density of states (DOS) at o = ath for the PM states of GdMg, GdZn and 
GdGd resolved into Gd (black curve) and Mg,Zn or Gd anion (red) components. The Gd d- 
component (green curve shaded underneath) and anion d-component (blue curve) are also shown. 
The upper (lower) panel shows the DOS for an electron spin-polarised parallel(anti-parallel) to the 
local moment on the Gd site. The total DOS, an average over all directions, is unpolarised. 


of the f-electron local moments as well as their effect upon it. Fig. shows the non-free 
electron-like PM valence density of states (DOS) of GdMg, GdZn and GdGd at ath for an 
electron spin-polarised parallel and anti-parallel to the local moment on the Gd site [TB] . 
Averaged over equally weighted moment orientations the DOS is unpolarised overall. Below 
Tc the electronic structure adjusts and spin-polarises [H] when magnetic order develops. The 
Gd f-moment interactions are properties of the electronic structure around the Fermi energy, 
ep. The Fermi surface (FS) of PM GdMg (for a = 0.96ath), Fig. |^a), shows a distinctive box 
structure so that a wave-vector, (0,0, ^), connects (nests) [MH3H] large portions of parallel 
FS sheets and drives AFl magnetic correlations. This topological feature is absent in GdZn’s 
and GdGd’s FS’s. Weak hybridization between Gd-5d and lower lying, nominally hlled Zn- 
3d or Gd-4d states, shown in Fig. causes complex differences between their electronic 
structures around ep and GdMg’s. In GdZn the Zn 3d bands are narrower than GdGd’s 
4d ones and lie at slightly higher energies [18]. Moreover we hnd that lattice compression 
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FIG. 4: The 3D Fermi surface for (a) PM GdMg (a = 0.96ath, where Fig.2(a) shows AFl order) 
which shows nesting, and (b) PM GdCd, (o = O.dAath, close to Fig.l(b)’s peak position) showing 
the ’hot spot’ at wavevectors k = (^, ^,0). The finite width of the FS features reflect the local 
moment disorder. 


increases Gd d-state occupation relative to sp-ones in these compounds [Ml ES] which affects 
FS topology. In particular, as shown in Fig. |^b) for GdCd, we find that Fig. [^b)’s peak 
correlates with a distinct electronic topological transition - a ’hot spot’ formed by a hole 
pocket around k = (|, |, 0), collapsing as a is reduced. 

Atomically localised f-electrons and their intricate physics is inevitably the focus for 
lanthanide material studies. But the valence electron glue in which the f-moments sit also 
harbors surprises. Its s-, p- and d-electrons can shift it far from a nearly free electron model, 
as exemplified by the canted magnetism of GdMg and the stark contrast of the magnetism 
of isoelectronic GdZn and GdGd with their disparate pressure variations. The predictive 
ab-initio computational modelling described here has successfully accounted for the subtle 
aspects of the valence electrons’ spin polarisability around Sp and how it is affected by 
occupation of lower-lying lanthanide-other metal d-electron bonding states. This implies 
that further successful quantitative modelling of the rich variety of technologically useful 
lanthanide-transition metal materials must also treat valence electronic structure accurately 
and in quantitative detail. We have shown that coordinated ab-initio theory-experimental 
studies have the capability of producing new guidelines for understanding the magnetism 
in lanthanide-transition metal magnets. Factors such as the average number of valence 
electrons or band-hlling, separation in energy of the lanthanide 5d and the other constituents’ 
d-bands and the valence band widths, reminiscent of the modern analogs of the famous 
Hume-Rothery rules HU for alloy phase stability, will influence the nature of the valence 
electron structure around ep and the magnetism it supports. 
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